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Abstract 

The electronic structure of the zig-zag bilayer strip is analyzed. The electronic spectra of the 
bilayer strip is computed. The dependence of the edge state band flatness on the bilayer width is 
found. The density of states at the Fermi level is analytically computed. It is shown that it has 
the singularity which depends on the width of the bilayer strip. There is also asymmetry in the 
density of states below and above the Fermi energy. 
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I. INTRODUCTION 

Carbon atoms can create a variety of forms such a graphite, diamond, carbon fibers, 
fullerenes and carbon nanotubes. A carbon nanotube can be described as a graphene sheet 
rolled into a cylindrical shape so that the structure is one-dimensional with axial symmetry 
and in general exhibiting a spiral conformation called chirality. They are interesting because 
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of their unique mechanical and electronic properties [lj. From the pioneering works [2|- 
the electronic properties of graphite have attracted interest because of unconventional phys- 
ical properties of a graphite layer. The development in the fabrication of the single layers of 
graphite (graphene) 5] caused a striking level of interest in the investigation of the carbon 
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compositions. In addition to the closed carbon molecules [6], systems with boundaries also 
show interesting features. The nanographite zig-zag ribbon possesses localized edge states 
near the Fermi level. States like that are absent for ribbons with armchair edges [7j. The 
graphite sheet is a zero-gap semiconductor with the density of states(DOS) vanishing at 
the Fermi level, the edge states of the zig-zag ribbons produce a peak in the DOS at the 
Fermi level. Both the carbon nanotubes and graphite layers have the edge states because 
of their boundary. The presence of the edge state results in the relatively important con- 
tribution to the density of states (DOS) near the Fermi energy. It was found [8|, |9[ that 
the HOMO-LUMO (highest occupied molecular orbital and lowest unoccupied molecular 
orbital, respectively) gap is inversely proportional to the length of the zig-zagcarbon nan- 
otube segment. The zig-zag ribbons have partly flat bands at the Fermi level [7|, llO{|. In the 
presented paper, we focus on the computation of the electronic spectra of the zig-zag bilayer 
strip and also on the computation of the DOS of the edge states near the Fermi level. 

II. THEORY 

Firstly, we describe the model for the zig-zag bilayer strip. We will study the edge and 
size effects using the tight-binding model for this strip shown in Fig. [TJ 

The 7r electronic structures are calculated from the tight-binding Hamiltonian 

i i,j i i,j 

+ J2 W m(\ ( Pt)(Vn\ + h.c), (1) 




FIG. 1: Structure of finite-width graphene bilayer. 

6i and el are the site energies of the upside and down layer; Iff), \<ff) are n the orbitals on 
site i at the upside and down layer; 7^, 7^ are the intralayer hopping integrals; Wy are the 
interlayer hoping integrals which depend on the distance dij and angle % between the 7Tj 
and ttj orbitals. 

To describe the parameter which characterizes the zig-zag bilayer strip, we start from 
the graphene layer where we can define the vectors connecting the nearest neighbor carbon 
atoms in the form: 
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r x = o(0; 



T2 = a{ 
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The distance between atoms in the unit cell is d = | 7$ | = 
to the double-layer graphene strip in the form: 

where 

M+l 
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-^. We want to find a solution 



^ u c?) = E (^^ + <^ Bi ; 



(3) 



(4) 
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and 
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¥{~?) = E ( C *M + C ^b) ■ (5) 

i=0 

We want to find a solution to the above equation in the form of the Bloch function 

Mt, ?) = -J= E e 4t -*W - ?„)>, (6) 



n 

where a denotes A or 5 atoms. Here fe = (/c x , 0), r n is the position of a unit cell and N is 
the number of a unit cell; \<f(r}) is a 7r orbital. We denote 

e t = e t = A = fo,«"(r - A,)| W'O" - ^)) = (<^(r - i^l^l^ - Sj)), (7) 
e 4 = ?i = -A = (v9 out (r - BOI^I^Cr - S*)) = ( V ou \r - A\)\H\y ou \r - A',)). (8) 

Now we define the intratube hopping integrals as jij = 7^ = 70. We take into account 
only the interaction between nearest-neighbors also in the case of interlayer interaction and 
denote 

( V {r - AftHMr - B'i)) = >n. (9) 

In confining the structure along the width, the edge states are induced by terminating the 
width dimension with zig-zag shaped edges. The presence of edges in the bilayer strip 
changes the dimensionality of the system from a two-dimensional to a one- dimensional sys- 
tem. The electronic spectrum of the zig-zag bilayer strip can be described by the following 
system of equations: 

(E - A)C Am = -7oC Bm _ 1 - 9kC Bm ~ 7iCj%., (10) 

(E + A)C Bm = - Jo C Am+1 - 9kC Am) (11) 

(E + A)C A>m = -70CV, - 9kC B > m , (12) 

(E - A)C B , m = -ioC A > m+1 - g k C A , m - 7l C Am , (13) 

where 

g k = 2j cos(k x a/2). (14) 

Here m = 1, , M, are site indices, where M describes the width of the graphene bilayer. 

We assume that the A and B M+1 sites are missing. So we have the boundary condition 
C Ao = C Bm+1 = C A i = C B > = 0. The solution is assumed to be [ill] 

C Am = Ae ipm + Be~ ipm , (15) 



C Bm = Ce ipm + De~ ipm . (16) 

Here A,B,C and D are the coefficients which have to be determined, and p is the transverse 
wave number. From the boundary condition we have 

C Ao = A + B = } (17) 

C Bm+1 = Ce ip{M+l) + De- ip{M+l) = 0. (18) 

And so 

C Am =A(e ipm -e~ ipm ), (19) 

C Bm = C(e ipm -z 2 e- ipm ), (20) 

where z = e tp ( M+1 \ And similarly, 

C A , m = A\e ipm -e- ipm ), (21) 

C B , m = C'{e ipm - z 2 e- ipm ). (22) 

Substituting Eqs.(fT9li2^1) into Eqs.(lTQll3"l) we obtain 

{E - A) (e iprn - z 2 e~ ipm ) A + [ 7o ( e ip( - m ~ 1) - e - ip{m - l) ) + g k (e ipm - e~ ipm )] C+ (23) 

+7l ( e *pm _ e -i P m^ C i = Q) 

[7o (e ip(m+1) - z V ip(m+1) ) + # fc (e ipm - z 2 e~ ipm )] A+(E + A) (e ipm - e~ ipm ) C = 0, (24) 

(£ - A) (e ipm - e~ ipm ) C + [ 7o (e^ (m - 1} - z 2 e -^ m ~^) + ^ ( e ^ m - z 2 e- ipm )} A'+ (25) 

+ 7l (e ipm - z 2 e~ ipm ) A = 0, 

[70 (e ip(m+1) - e" ip(m+1) ) + g k (e ipm - e~ ipm )] C + (E + A) (e ipm - z 2 e - ipm ) A' = 0. (26) 

This homogenous system of equations has a solution only if the following conditions are 
fulfilled: 

[E 2 - A 2 - 7l (£ + A) - ( 7o e- jp + g k ) ( lo e vp + g k )] e 2%pm + 

z * [E 2 - A 2 - 7l (£ + A) - ( 7o e-* p + g k ) ( 7o e ip + g k )] e~ 2ipm - 
- (E 2 -A 2 - lx {E + A)){z 2 + 1) + (g n + 7o e ip ) 2 + z 2 (g n + lo e~ ip ) 2 = 0, (27) 

or 

[E 2 - A 2 + lx (E + A) - ( 7o e- ip + g k ) ( lo e ip + g k )] e 2ipm + 
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z 2 [E 2 - A 2 + 7l (E + A) - ( 7o e-* + g k ) ( 7o e* p + &)] e" 2 ^- 

- (£ 2 - A 2 + 7l (£ + A))(^ 2 + 1) + (g k + 7o e^) 2 + ^ 2 (<? fe + 7o e"^) 2 = 0. (28) 

The coefficient of e ±2pm terms and the constant term have to be equal to zero. Thus, we 
obtain the energy spectrum 



E h2 = ^ ± J 7o 2 + 2 im cos(p) + g\ + (^ + A) , (29) 



£ 3 ,4 = -y ± yj 7o 2 + 2 7 o^ cos(p) + <£ + (^ - a) . (30) 

and the equation which gives the longitudinal wave number p is 

sin \pM] + — sin \p(M + 1)] = 0. (31) 

7o 

For M 3> 1 Eq. (l3TI) can be written as 

sin [pM] = 0. (32) 

The solution is given by 

P=|/. (33) 

Substituting this solution into Eq.f l29|30|) we get the energy spectrum for the graphene 
bilayer with the periodic boundary condition along the y-axis. 

III. EDGE STATES OF GRAPHENE BILAYER 

Now we are interested in the edge state of the graphene bilayer. This solution can be 
obtained in the form p = tc + ir/ 12J. We get the following equation for rj: 

sinh [r]M] - — sinh [r](M + 1)] = 0. (34) 

7o 

The edge state can exist when the condition 

\2cos(k w a/2)\< 1+ 1 : (35) 

is fulfilled. The energy spectrum of a state like that is given as 



E 1>2 = ^ ± J 7o 2 - 2 lo9k coshfa) + g 2 + (^ + a) , (36) 



£s,4 = -y ± Jl 2 - 2 l0 9k coshfa) + gl + (^ - a) . (37) 



For big enough M the solution of Eg. (134]) can be expressed in the form 12 ] 

7] = In 



i-4 

°k + 2M+1 



where 1/ck = |2cos (^f) |. We denote k x = k. From Eq.( l38|) we have 

1 + cl (c|_- 1 " 
\ 
and so 



(38) 



cosl1 ^^- ^ > (39) 



^ = #±A7g^S^+ ^+A , (40) 



H 


/TO c 2M+4 + U J ' 


-y ± 





^=-f ± v«>ri L+ if- A ;- (4l) 

Now we assume, similarly as in [13|, that 71 > 2 A and also it is assumed that the width 
of the graphene bilayer is big enough and the following condition is fulfilled: 

(c 2 - 1? 

7i » 7o 2 Hto L - ( 42 ) 

c k 

The bands are given by 



,,2 („2 -\ 2 



£ lW ^-» + A + 7 J° 2A ( ^ M ; 4 ' , (43) 

^^- Ay <44) 

^*)=- a+ a^' <45) 

B 4 (*) = -71 + A - j^T-%^ - (46) 

We are interested in the E 2 (E 3 ) band which is the valence (conductance) band of the edge 
states. The minimum of the E 2 band is 

E2 > mm = ~ A ~ ^M 77T2A' (47) 



where it was used that e x = (1 + x/M) M for M — > oo and also Eq. (l35|) . We found E 2>m ax = 
—A. Similarly for the E% band 

£ 3 , ma , = -A + e2M ^_ 2A , (48) 

and E^^min = — A. We can see that the width of the band is inversely proportional to the 
width of the bilayer. 

The density of states can be expressed in the form 

N(E) = ~, (49) 

dk 

where L is the length of the bilayer in the x direction. We get the density of the state in 
the region E^^ m i n — E^ max in the vicinity of the energy E = —A in the form 

N W = ^^ iT^TTTT- ( 5 °) 

2vraM (E + A)^+^ ( 7o 2 /( 7l - 2A)) 1/2(M+1) 

The density of the state in the region E 2jm i n — E 2tma x i n the vicinity of the energy E = — A 

has the form 

N{E) = 2naM (E + A)^ ( 7o 2 /(7i + 2A)) 1/2(M+1) ' (51) 

Both these densities of the states have a singularity at the energy E = —A. The strength 

of the singularity depends also on the width of the bilayer. The width of the bilayer is 

characterized by the parameter M. The density of the state of the E4,(k)(Ei{k)) band is the 

same as the E 3 (k)(E 2 (k)) band. 

IV. CONCLUSION 

In the presented paper the electronic spectra of the zig-zag bilayer strip was studied 
analytically. We get for big enough M that the electronic spectra of the graphene bilayer 
strip are similar to the spectra of the graphene bilayer with the periodic boundary condition. 
Because of the boundary we also get edge states. It was shown that the width of the edge 
state band is inversely proportional to the width of the bilayer strip which is characterized 
by the parameter M. So for big enough M we get partly flat bands of the edge states. The 
density of states at the Fermi level has a singularity which also depends on the width of the 
bilayer strip. There is asymmetry in the DOS at the Fermi energy, similarly to the electron- 
hole asymmetry in the bilayer graphene [l4j. This asymmetry is caused by the parameter 
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A which describes the difference in the site energy of the atoms at the sites Ai, B i and the 

atoms at the sites Bi,A { . 
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